Electron as soliton: Nonlinear theory of dielectric polarization 
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We describe a class of theories of dielectric polarization and a species of solitons in these theories. 
The solitons, made entirely out of the polarization field, have quantized values of the electric charge 
and can be interpreted as electrons and holes. A soliton-antisoliton pair is an exciton. We present 
numerical evidence that the elementary soliton is stable. 



Absorption of light by dielectrics is one of the oldest 
problems in quantum theory. Because a typical frequency 
of visible light is much larger than the maximal frequency 
of phonons, absorption of even one photon must result in 
many phonons. Frenkel's idea [l[ was that the immense 
suppression, naively expected for such a process, can be 
overcome if there is an intermediary — an electronic ex- 
citation of the solid (an exciton). An alternative to this 
point of view would be to imagine that phonons pro- 
duced by light absorption initially form a highly coherent, 
nearly classical state. Then, the effective nonlinearity is 
enhanced by the large amplitude of the phonon field in 
the final state, resulting in a sizable absorption rate. 

In the present Letter, we would like to argue that these 
two points of view are, in fact, complementary to each 
other. Namely, if the phonon field is treated as a nonlin- 
ear system, there are solitons of this field that are nat- 
urally interpreted as electrons and holes, and a soliton- 
antisoliton pair can be interpreted as an exciton. 

We wish to stress from the outset that our results are 
quite different from the familiar polaron picture. In the 
latter case, there are two fields, which describe electrons 
and phonons separately, and the polaron is a bound state 
involving both. In our case, there is no separate elec- 
tron field — electrons emerge as nonlinear excitations of 
phonons, somewhat similarly to how a baryon emerges 
as a nonlinear excitation of pions in Skyrme's model 
0. An even closer similarity is with the solitons of the 
sine-Gordon (SG) model, which can be viewed [3, 0| as 
fermions of the massive Thirring model. 

Neither of these similarities is complete, however. Both 
in the SG case and in Skyrme's model, solitons exists 
directly in the continuum limit, while in our case, at least 
in the present version of the theory, they can only live on 
a lattice. 

One reason to pursue the picture of electrons as soli- 
tons is that it allows one to think about various processes 
of interest in nanoscale dielectrics entirely in classical 
terms. For example, one may want to know if a highly 
excited excitonic state relaxes to the lowest-energy ex- 
citon via carrier multiplication (i.e., production of more 
excitons) or via emission of multiple phonons. By view- 
ing electrons as solitons of the phonon field, we obtain a 
classical approximation, which can be used to simulate 
this process on a computer. 

Quantum physics of these solitons, in particular, the 



question of their spin and statistics, also promises to be 
interesting. In the present Letter, however, we limit our- 
selves to problems of the classical theory: arguing the 
existence of solitons and describing their simplest prop- 
erties. While a plausible existence argument can be given 
analytically, as further evidence we also present results 
of a numerical study. 

Consider then the field 0„(x, t) of an optical phonon. 
We will consider in parallel spatial dimensions d — 2 and 
3 (planar and bulk dielectrics, respectively). The number 
of components of 0„, i.e., the number of values that n 
takes, is equal to d. We assume that 6'„ is defined on 
the faces (for d = 3) or edges (for d = 2) of a rectangular 
lattice with spacings a„. Up to a constant factor. On is the 
component of dielectric polarization in the n-th direction. 
For notational simplicity, we will often use continuum 
notation for lattice derivatives. Thus, for example, the 
electric charge density (in suitable units) will be written 
as a divergence, 

p(x,t) = V-6>(x,t) = a„0„(x,i) , (1) 

but this should be read to mean 
d ^ 

p(5c,i) = V— [0„(x + n,i)-0n(x,<)] , (2) 

where n is the primitive vector in the n-th direction. Note 
that, since 0„ live on the faces (or edges), p naturally lives 
in lattice cells (i.e., on the lattice dual to the original). 
Hence the tilde in eq. 

The form of the Lagrangian density that governs dy- 
namics of 9n should ultimately be based on the properties 
of the material. For our present purposes, it will be suf- 
ficient to use the simplest form, which we write as the 
first few terms of a derivative expansion: 
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L = -v{e)--{w-ef 
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(3) 

where C > 0, £, and Mmn are constants. In what fol- 
lows, we limit ourselves to searching for static solutions 
of the theory ([S]), and for these the terms containing time 
derivatives do not matter. They would be crucial, how- 
ever, if we were to study the phonon spectrum. 

The potential V contains no derivatives of 0„ and de- 
scribes the dielectric response of the medium to uniform 
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static fields. The only feature of it essential here is 
a certain periodicity, so we use the simple cosine form 
(which in one dimension would be the potential of the 
SG model): 



(4) 



In ([4]), /i'^ is a parameter, and (t„ is the area of the faces 
(for d = 2, length of the edges) that are orthogonal to 
the n-th direction. At small On, V « ^is^9^, which allows 
one to relate /i^ to the bulk dielectric constant e: fi^ — 
e^/7r(e — 1), where = 14.4 eV A is the electron charge 
squared. 

The second and third terms on the right-hand side of 
([3|) describe the capacitive and inductive effects due to 
short-scale charge separation. In particular, 



Ip = 1/fiVC 



(5) 



is a characteristic length for spatial variations of the po- 
larization charge. For a material with a high degree of 
spatial symmetry, where we expect the expression ([3]) to 
apply, C, together with C and Mmm can in principle be 
determined by fitting the phonon spectrum. 

Note that if /i is small (i.e., e is large), as in many 
semiconductors, Ip can be much larger than the atomic 
scale. In this case, in nanocrystals of sufhciently small 
size, effects of spatial dispersion may be quite significant. 
On the experimental side, deviations of measured exci- 
tonic spectra from predictions of the "envelope" theory 
(in which electrons are considered point-like) have been 
reported for PbSe crystals with diameters smaller than 7 
nm Q. 

The short-scale effects encoded in the values of C, £, 
and Aimn should be distinguished from effects due to 
macroscopic electromagnetic (e.-m.) fields that can exist 
in a dielectric. To take those fields into account, we would 
need to add to ^ the term 
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AqS/ ■e + -A-dte 



(6) 



where Aq and A are the scalar and vector potentials, 
e is the electron charge, and c is the speed of light in 
vacuum. In d = 3, this term also contributes to the 
phonon frequency. For the moment, we are not interested 
in macroscopic e.-m. fields of solitons, so in what follows 
we omit the Lagrangian ([6]) . Effects of impurities are also 
readily incorporated into the theory, by adding a term 
ti(x)V • 0(x, t), where w(x) is the impurity potential. For 
now, these will be neglected as well. 

An important property of eq. ([3]) is the absence of 
"transverse capacitances" , i.e., terms of the form {dmOnY 
with m ^ n. As we will see, such terms lead to a linear 
(confining) potential between solitons. Since our goal 
is to describe a material in which electrons and holes 
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FIG. 1: The original lattice (thin lines) and closed and open 
strings on the dual lattice. 



can be separated away to large distances, these terms 
are in fact forbidden. On the other hand, we could add 
terms of the form [dm sin(cr„0„)]^. At small 0„, these are 
indistinguishable from the forbidden terms, but at larger 
On they display the periodicity property that causes the 
tension of the "string" connecting solitons to vanish. If 
these terms are relatively small, then adding them to ([3]) 
will only deform our solution a little, and for the present 
work we have left them out. In applications to specific 
materials, they can be readily included. 

The argument forbidding the "transverse capacitan- 
ces" can be made more formal by identifying a relevant 
symmetry. For definiteness, we will speak about the 
d = 2 (planar) case; the argument for d = 3 is entirely 
similar. Consider a directed closed contour on the dual 
lattice (see Fig. [1]). Imagine changing the values of On 
at all the edges (in d — 3, faces) the contour crosses 
by ±27r/tT„, the sign depending on whether the contour 
runs parallel or opposite to the n-th axis. This transfor- 
mation will be referred to as adding a closed 27r string. 
The symmetry in question is the requirement that under 
this transformation the energy does not charge. 

An open 2tt string, such as the one shown in the lower 
portion of the figure, does cost energy, i.e., it is an ex- 
citation. If the symmetry with respect to adding closed 
strings is in place, however, the string tension is zero, and 
the energy of an open string is accumulated only near the 
ends. Since a 2tt string carries electric flux of 2tt, the ends 
of an open string are charged, with quantized charges 
equal to ±27r, so this excitation can be interpreted as an 
electron- hole pair. 

In the one-dimensional SG model, this way of creat- 
ing excitations leads to a soliton-antisoliton pair that is 
nearly stable (the more so the larger is the distance be- 
tween the solitons). Indeed, imagine starting from the 
state with On — (in one dimension. On has only one 
component) and then changing On by 27r at all sites of 
a chosen segment. If one then "cools" down the system, 
to reach a low-energy state, the field will become smooth 
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at the ends of the segment but will remain equal to 2tt 
in the middle: due to the high potential barrier, it can- 
not unwind back to zero there. Similarly, in our case, 
we also expect formation of a stable soliton pair (and, as 
we discuss shortly, this has been confirmed in numerical 
simulations). 

If, however, we have allowed the "transverse capac- 
itances", the symmetry with respect to adding closed 
27r strings would be broken. For an open string, energy 
would then have come not only from the ends but from 
the entire length and, for a large length, would be propor- 
tional to it. As a result, a pair of well separated soliton 
and antisoliton would not be an approximate solution: 
the linear potential would pull the soliton and antisoli- 
ton towards each other until they annihilate. (This has 
also been confirmed numerically.) 

The symmetry with respect to adding closed 2t: strings 
is local, in the sense that the number of independent 
transformations scales linearly with the volume of the 
system. Indeed, any closed loop can be viewed as a su- 
perposition of elementary closed loops, each encircling 
just a single site (in d = 2) or edge (in d = 3) of the 
original lattice. With the loop's directionality taken into 
account, there are twice as many ways to add an elemen- 
tary loop as there are sites or edges, respectively. 

The local character of the symmetry suggests that the 
configuration space of our theory is not the space of 
all possible 0„(x), but the space of equivalence classes 
with respect to adding an arbitrary number of closed 
27r strings, i.e., states that differ only by the number 
of closed 2tt strings should be considered as one state. 
The role of this equivalence relation will become clearer 
shortly, when we consider the limit in which our soliton 
turns into a vortex or a magnetic monopole. However, 
this relation is not essential in the method we have used 
for finding solitons. 

The presence of a string, similar to Dirac's string, 
suggests that our solution is related to the magnetic 
monopole (or the vortex in d = 2). Indeed, it is a gener- 
alization of those. The monopole is obtained in the limit 
C — *■ 0. Let us first consider this as a formal limit, i.e., 
as the requirement that V • = 0. Then, 9 is the curl 
of some other field, (/). In d = 2, this field lives on the 
sites of the original lattice and has only one component, 
so that 
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(e„m is the antisymmetric unit tensor); in d 3, lives 
on the edges. Changing (f> by 27r at a single site creates, 
in terms oi 6, an elementary closed 27r string encircling 
that site. So, imposing the equivalence relation as above 
amounts to requiring that cj) is an angular variable with 
period 27r. Substituting eq. ^ into eq. ([4]), we see that 
the potential V{d) becomes the energy density of the 
2-dimensional XY model, which has well-known static 
solutions — vortices. In d = 3, a similar argument leads to 



compact electrodynamics, whose solutions are magnetic 
monopoles @ . (They are "magnetic" with respect to the 
gauge field that replaces the scalar in the curl condition, 
not with respect to ordinary electromagnetism.) 

A less formal way to view the C — > limit is to note 
that in this case the screening length Ip, eq. (O, goes to 
infinity. Since the total charge of an elementary soliton 
is fixed (and equal to 2tt), the charge density V • 6 goes 
to zero. Thus, if a soliton and an antisoliton are sepa- 
rated by a distance much smaller than Ip, i.e., sit well 
within each other's polarization cloud, we expect them 
to interact as if they were monopoles (or vortices) of the 
respective C = theories. 

As well known, monopoles and vortices have long- 
range interactions, which are important in various prob- 
lems of statistical physics From our present stand- 
point, these interactions are a feature of the Ip —^ oo 
limit. For a finite Ip, the interaction range is of order Ip 
(provided, of course, that the c.-m. interactions ((6|) have 
been switched off). 

We now turn to numerical results. To search for static 
solutions of the theory ([3]), we drop the terms with time 
derivatives and numerically solve the relaxation equation: 



de 



1 dV 



Pp do 



(8) 



where r is a fictitious time. To avoid influence of bound- 
aries, we have used periodic boundary conditions. 

As the initial condition, we have used an open string, 
such as shown in fig. I, with various values of the elec- 
tric flux. The actual solutions (endpoints of the relax- 
ation process) are expected to have the flux quantized 
in integral multiples of 2-k — these are the values of fT„6'„ 
at which the potential ^ has minima. It is interesting, 
however, to start with an arbitrary value and see how 
the nearest quantized value is approached. Indeed, some 
of the most illustrative results are obtained by starting 
with UnQn close to a maximum of the potential. 

Fig. 2 shows results from a 2-dimensional 64x 32 lattice 
for the case when the initial string has a flux of ci^i = 
3.027r; results from a 3-dimensional 64 x 32^ lattice (for 
the same initial flux) are similar. The screening length 
{Ip = 1) in either case was about 3 times the lattice 
spacing (which was the same in all directions). 

The initial flux is just larger than 37r, a value at which 
the potential has a maximum. So, we expect the flux 
(and therefore also the charge at the string's end) to roll 
to the nearest stable value, aiOi — An, which is twice the 
minimal soliton charge. The final value of the charge at 
each end is indeed close to ±47r (a small difference from 
±47r being attributable to finite-size effects), but we see 
that it is achieved by forming two solitons at each end, 
rather than a single doubly-charged one. The final state 
was stable on the timescale of our simulation. 

We have observed formation of two solitons instead 
of one even when the initial string flux already was 
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FIG. 2: Surface plots of the charge density towards the be- 
ginning (top) and end (bottom) of the relaxation process. 



(7iOi = Att. Wc interpret this as an indication that the 
short-range interaction between sohtons of like charge is 
repulsive, and a compound soliton is unstable with re- 
spect to decay into elementary ones. 

To summarize, we have described a class of theories 
that provide a unified description of phonons and exci- 
tons in a dielectric. Phonons are small fluctuations of 
a polarization field, and excitons are soliton-antisoliton 
pairs made of that same field. We have presented numer- 
ical evidence that the elementary soliton is stable and 
identified the symmetry (adding closed 2ii strings) that 
underlies both this stability and the quantization of the 
soliton charge. 
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